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Abstract. Let a,b £ Q* be rational numbers that are multiplicatively independent. 
(~^\ . We study the natural density <5(a, b) of the set of primes p for which the subgroup 

of F* generated by (a mod p) contains (b mod p). It is shown that, under assumption 
of the generalized Riemann hypothesis, the density 5{a, b) exists and equals a positive 



rational multiple of the universal constant S = Y[ p pr i me (l ~p/(p 3 ~ !))• An explicit 
value of <5 (a, b) is given under mild conditions on a and b. This extends and corrects 
earlier work of Stephens [13]. 
■ We also discuss the relevance of the result in the context of second order linear 



recurrent sequences and some numerical aspects of the determination of S(a,b). 



1. Introduction 

O 

Artin's original conjecture on primitive roots gives, for each non-zero integer a, a 
C\) \ conjectural value 5(a) of the density of the set 

On 

Q\ • (1-1) {p prime : (a mod p) = F*} 

^ - 

inside the set of all primes. It equals 5(a) = c a ■ Y\ p prime (l — p ( p -i) )) where c a is 
some explicit rational number that is positive whenever a is not equal to —1 or to 
a square. Artin's conjecture was proved by Hooley [3] under the assumption of the 
generalized Riemann hypothesis. Unconditionally, there is not a single value of a 
for which the set of primes in (1.1) has been proved to be infinite (cf. [9]). 
In this paper we study the density 5(a, b) of the similarly defined set 

(1.2) {p prime : b mod p G (a mod p) C F*} 

inside the set of all primes. In view of our application in Section 6, we allow a and b 
to be non-zero rational numbers, and exclude the finitely many primes dividing the 
numerators and denominators of a and b from consideration in (1.2). 

If a and b satisfy a multiplicative relation a x b y = 1 for exponents x, y G Z that 
are not both equal to zero, then one can prove unconditionally that 5(a, b) is a 
rational number, and that it is positive in all but a few trivial cases [14]. We will 
therefore restrict to the case where a and b are multiplicatively independent in Q*, 
i.e., no non-trivial relation of the type above holds. In this case, the following 'two 
variable Artin conjecture' has been proved unconditionally. 
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Theorem 1. Let a, b G Q* be multiplicatively independent. Then the set of primes 
defined by (1.2) is infinite. 

Theorem 1 is actually a special case of a theorem of Polya [10], and we include its 
short and elementary proof in Section 6. It does not show that the set (1.2) contains 
a subset of primes of positive density. 

We will mainly be concerned with the density of the set (1.2). Here the basic 
result is the following. 

Theorem 2. Let a, b G Q* be multiplicatively independent, and assume the validity 
of the generalized Riemann hypothesis. Then S(a,b) exists and equals 

S(a,b) = c a , b . J] a-^i) 

p prime 

for some positive rational constant c a ^. 

The constant c a j> in the theorem depends on the degrees of the number fields 
(1.3) i ! ij = Q(Cy,a 1/y ,6 1/< ) 

for i,j G Z>o. Here Qj denotes a primitive ij-th root of unity. As an explicit formula 
for c a ^ is rather cumbersome to write down, we only compute it explicitly in the 
'generic case' where the factor group Q*/ (—1, a, b) is torsionfree. This means that 
±a x b y is not an n-th power in Q* when x and y are not both divisible by n. 

Theorem 3. Let a, b G Q* be multiplicatively independent, and suppose that the 
group Q*/(— l,a, 6) is torsionfree. Define r(n) for n G Z^ by 

. _r,4-3 ordp(n) 

V 1 , 11 p 3 -p- 1 

p\n prime 

Then the constant c a ^ in Theorem 2 has the value 

c a>6 = 1 + r(lcm(2, A(a))) + e(6)r(A(6)) + e(ab)r(A(ab)). 

Here A(x) G Z denotes the discriminant of the quadratic field Q(y / x) for x G Q* ; 
and we put 

e(a . ) = |^ */A(;r)z 50 dd; 

! 1 if A(x) is even. 

The universal constant S = Yl prime (l — p/(p s — 1)) in Theorem 2, which is the 
analogue of Artin's constant A = Yl pr i me (l — Vp(p — 1)) arising for the original 
Artin conjecture, already occurs in Stephens' paper [13]. Our Theorem 2 occurs for 
positive coprime integers a and b that are not perfect powers as [13, Theorem 3], but 
the explicit value for c a ^ given there is involved and incorrect. The analytic part of 
Stephens' proof, which we summarize in the next section, is correct and generalizes 
in a rather straightforward way to our more general situation. His explicit evaluation 
of c ay b however, which is only carried through in one out of the eight subcases 
distinguished in [13], is incorrect, yielding an expression that is symmetric in a 
and b. Our proof of Theorem 3 separates the elementary calculus of double sums 
from the algebraic facts concerning the field degrees [F it j : Q]. 

In the final two sections we address the relevance of our results in the setting 
of second order recurrent sequences and deal with some numerical aspects of the 
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2. Results of Hooley and Stephens 

The proof of the special case of Theorem 2 occurring in [13] proceeds along the 
lines of Hooley's proof [3] of the original Artin conjecture. 

Artin's basic observation is that, for a G Q* arbitrary and p a prime number with 
ordp(a) = 0, the index [F* : (a)] is divisible by j if and only if p splits completely 
in the splitting field Fj = Q(C J ,a 1 / J ) of the polynomial — a over Q. By the 
Chebotarev density theorem, the set of these primes has natural density l/[Fj : Q]. 
The primes p for which a is a primitive root modulo p are those primes that do not 
split completely in any of the fields Fj with j > 1. In fact, it suffices to require that 
p does not split completely in any field Fj with j prime. For Artin's conjecture, a 
standard inclusion-exclusion argument readily yields the heuristic value 

The right hand side of (2.1) converges whenever a is different from ±1, since in that 
case [Fj : Q] differs from its 'approximate value' j ■ (p(J), with ip denoting Euler's 
(^-function, by a factor that is easily bounded in terms of a. In fact, we obtain an 
upper bound for the upper density of the set (1.1) in this way. In order to turn this 
heuristic argument into a proof, Hooley employs estimates for the remainder term 
in the prime number theorem for the fields Fj that are currently only known to 
hold under assumption of the generalized Riemann hypothesis. 

In the situation of Theorem 2, one can find a conjectural value for S(a,b) in a 
similar way. For each integer % > 1, one considers the set of primes p with ord p (a) = 
ord p (6) = for which the index [F* : (a)] is equal to i and the index [F* : (b)] 
is divisible by i. These are the primes that split completely in the field F^i = 
Q(C», a 1 / 1 , b 1 / 1 ), but not in any of the fields F id = Q(Cij, a 1/iJ , b 1 / 1 ) with j > 1. As 
before, inclusion-exclusion yields a conjectural value for the density Si(a,b) of this 
set of primes, and summing over i we get 

OO OO OO , .n 

(2-2) *M) = £^M)££r/^- 

i=l i=l j=l I- *'J • ^ 

Note that Si(a,b) is nothing but the primitive root density 5(a) from (2.1). 

As in the case of (2.1), the right hand side of (2.2) converges if the degrees 
[Fij : Q] are not too far from their 'approximate values' i 2 j ■ (p(ij). As we will see 
in Lemma 3.2, this is exactly what the hypothesis that a and b be multiplicativily 
independent implies. 

The proof of (2.2) by Stephens [13], which assumes the Riemann hypothesis 
for each of the fields F it j, closely follows Hooley's argument in [3]. The restrictive 
hypotheses on integrality and coprimality of a and b are not in any way essential. 
The only requirement for the argument to work is that, up to a factor that can 
be uniformly bounded from below by some positive constant, [F i} j : Q] behaves 
as i 2 j ■ <p(ij). We will show this in the next section, so there is no need for us to 
elaborate any further on the proof of (2.2); we will merely be dealing with degrees 
of radical extensions in order to prove Theorems 2 and 3. 

The universal constant S = Yl prime (l — p/(p 3 — 1)) is the value of the right 
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constant A = f| pr i me (l ~~ ~~ 1)) is obtained from the right hand side of (2.1) 

by putting [Fj : Q] = j ■ <p(j). The 'correction factors' c a in Artin's conjecture and 
c a ,b in Theorem 2 measure the deviation of the field degrees [Fi : Q] and [Fij : Q] 
from these values. As in the case of Artin's conjecture, the basic problem is that 
radical extensions of Q involving square roots are not in general linearly disjoint 
from cyclotomic extensions. In the case of Theorem 3, the situation is sufficiently 
simple to allow an expression for c a ^ by a formula that can be fitted on a single 
line. As the proof of Theorem 2 will show, all other cases can in principle be dealt 
with in a similar way. 

For each prime q, the corresponding factor 1 — l/q(q — 1) in the definition of 
A has a well known interpretation: it is, generically, the fraction of the primes p 
for which (a mod p) generates a subgroup of F* of index not divisible by q. In a 
similar way, the factor 1 — q/(q 3 — 1) at a prime q in the product for S represents, 
generically, the fraction of the primes p for which the number of factors q in the 
index [F* : (a mod p)] is at most the number of factors q in [F* : (b mod p)]. 

As in the case of Artin's conjecture, we cannot prove unconditionally that the set 
of primes in (1.2) has positive lower density. The problem is that Chebotarev's den- 
sity theorem only allows us to simultaneously impose conditions of the type above 
at primes q for finitely many primes q. As the fraction of primes that is eliminated 
by imposing a condition at q is generically positive, one can prove unconditionally 
that there exists a set of primes of positive density for which b mod p is not con- 
tained in the subgroup (a mod p) of F*. This was already noted by Schinzel [12]. 
For the set (1.2) itself, the best unconditional result available is Theorem 1. 

3. Radical extensions 

For Artin's original conjecture, one has to compute the degree of Fj = Q(£j, a 1 ^), 
the number field obtained by adjoining all j-th roots of an element a G Q* \ {±1} 
to Q. The result may be found in [15, Prop; 4.1]. The key observation is that if we 
take a ^ ±1 such that Q*/ (—1, a) is torsionfree, then a is a square in Q(Cn) if an d 
only if the discriminant A(a) of Q(y/a) divides n; moreover, a is not a k-th power 
with k > 2 in any cyclotomic extension of Q. If j\ divides j and a is as above, then 
Kummer theory yields 

(3.1) [Q(0,a 1/J1 ):Q] = 

For arbitrary a G Q*\{±1} and j G Z>o, let t = gcd{ord p (a) : p prime} be the order 
of the torsion subgroup of Q*/(— l,a), and take j± = j/ gcd(j, t). There are two 

cases. If a is a t-th power in Q*, say a = a\, we have Fj = Q(Q, a 1 /- 7 ) = Q(Cj, o}^ 1 ) 
and (3.1) can be applied directly. If a is not a t-th power, then t is even and —a = a\ 
is a t-th power in Q*. We have an extension 

Q(Q, Gya!'*) = Fj C F> = Q(C 2j , a^' 1 ) 

of Fj of degree at most 2 for which the degree [Fj : Q] is given by (3.1), and one 
is left with the determination of [Fj : Fj] G {1,2} as in [15]. A somewhat subtle 
case distinction is necessitated by the peculiarity that the element —4, which is 



j iji ■ (p(j) if 2 divides ji and A(a) divides j; 
I ji ' pU) otherwise. 
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have [Fj : Fj] = 2, then j and t are even and 2ji divides j. Thus, the 'degree loss' 
j<p(j)/[Fj : Q] with respect to the generic value j<p(j) is always an integer dividing 
2t. The factor 2 reflects the fact the kernel of the natural map 

Q7Q* fc -> Q(C*)7Q(C*)** 

is an abelian group annihilated by 2. More precisely, it vanishes if k is odd; if k is 
even it is generated by the elements x k / 2 Q* satisfying A(x)\k and, for k = 4 mod 8, 
the element -2 fc / 2 Q* fc . 

In our two variable setting, where we deal with the fields Fij = Q(Cij, a 1//jJ , b 1 ^), 
the statement in terms of the map above easily leads to the following generalization. 

3.2. Proposition. Let a, b G Q* be multiplicatively independent, and let t be the 

order of the torsion subgroup of Q*/ (—1, a, b) . Then for all i,j G Z>o, the quantity 

JiJ 



is a positive integer dividing At. In the torsionfree case t = 1, it is equal to 

- the number of elements in {1, A(a), A(&), A(ab)} dividing ij if i is even; 

- the number of elements in {1, lcm(2, A(a))} dividing ij if i is odd. 

Proof. Pick i,j G Z>o, and write k = ij. Let Wij C Q* be the subgroup generated 
by a and V = b k l % ', and Wij the image of Wij in Q*/Q* fc . As the order of W itj 
divides ik = i 2 j, we write i 2 j = #Wij ■ tij with tij G Z>o- By Kummer theory, 
the degree of 

= Q(Cy»a 1/iJ '»6 lA ) = Q(Cfc, V 7 ^) 

over Q(Cfc) equals #V , [^ij']) with ip : Wij — > Q(Cfc)*/Q(Cfc)* fc the natural map. 
We deduce that the 'degree loss' foj for F^j can be written as foj = tij ■ # ker ^. It 
is a decomposition of /j^ into a factor £jj coming from 'torsion in Q*'and a factor 
#ker-0 measuring the additional torsion caused by the adjunction of (k- 

As Wij is a finite abelian group on 2 generators and ker^ C Wij is annihilated 
by 2, it is clear that ^kertfj divides 4. In order to show that tij divides t, we let 
T C Q* be the inverse image of the torsion subgroup of Q*/(— l,a,b) under the 
natural map Q* — > Q*/ ( — 1, a, 6). Then T contains F = (—1, a, 6) as a subgroup of 
index £, and W \j is the subgroup of T/T k C Q*/Q* fc generated by a and 6 J = 
The integer t^j is the order of the kernel of the composed map 

(^V) W *,3 ^vlV k ^TlT k 

{a k ,b k )WijDV k 1 1 ' 

so it divides #ker [V/V k -> T/T fc ] . As and T/T fc are finite abelian groups 

of the same order, #ker/ = #coker/ = [T : VT k ] divides [T : V] = t. This shows 
that fij divides At. 

Assume now that Q*/ (—1, a, b) is torsionfree. Then we have T = V and Uj = 1 
in the argument above, and fij = #ker^. Clearly, ker if; vanishes if k is odd. If k 
is even, the 2-torsion subgroup of Wij is generated by a fc / 2 if i is odd and by a k l 2 
and 6 fc / 2 if i is even. As ip vanishes on the residue class of x k l 2 G (a fc / 2 , 6 fc / 2 ) in Wij 
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3.3. Corollary. The integer fcj in 3.2 only depends on the greatest common di- 
visors gcd(z,2f) and gcd(ij,8st), where s is the product of the primes p for which 
ordp(a) and ord p (b) are not both equal to 0. 

Proof. In the proof above, one needs gcd(fc, t) = gcd(zj, t) to determine the kernel 
ker [V/"l/ fc — > T/T k ~\ and gcd(z,£) to determine the order tij of the intersection of 
Wij = (a, b k / 1 } with this kernel. As A(x) divides 4s for all x G V, we can determine 
the kernel of V/V k -> Q(Cfc) VQ(Cfc)* fc if we know gcd(fc, 8s) = gcd(zj, 8s). Knowl- 
edge of the parity of % enables us to intersect this kernel with W^j/iW^j fl V k ), 
thus yielding the second factor #kerijj in tij. □ 

4. Evaluation of the basic double sum 

From (2.2), 3.2 and 3.3 it is clear that, in order to evaluate 8(a, b), we need to 
evaluate for non-zero integers m, n the double sum 



oo oo 



(4-1) s m , n = J2H 



Mi) 

m\i mn\ij 



This is a rather straightforward computation in elementary number theory leading 
to the following result. 

4.2. Theorem. Form, n e Z >0 , the value S m>n of the series in (4.1) is the rational 
multiple 

g TT -P 4 TT P 3 +P 2 
Om.n o o I I o ills 1 

J- J- P J — p — 1 J- J- — p — 1 



of the universal constant S = Y[ p pr ime(l ~~ p^i) occurring in Theorem 2. 

Proof. We can sum over all i > 1 in (4.1) after substituting mi for i. Putting 
ij = nd and summing over all d > 1 then yields 

1 °° 1 

d=l ^ v y j|nd 

Writing J = n p |xP f° r ^ ne largest squarefree divisor of x, we have for any integer x 
$^J>(j) = ^3V(j)v(x) ^2jp(x/j)n(x)(p(x). 

j\x j\x j\x 

This enables us to write 

s 1 n(nd)<p(nd) fi(n)<p(n) 

m ' n m 2 n 2 ^— ^ d 2 (p(mnd) m 2 n 2 (p(mn) ^—J 

where / is the multiplicative function defined by 

,s _ 1 <p{mn) n(nd) <p(nd) 

J \ U J to " 7 tT " 7~\ " 7~A * 
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As YldLi /(^) ^ s absolutely convergent, we can use the values 

{_pi-3k for p { mn; 

p-3k £ Qr p| n . 

— (p — l)p~ 3k for p|m, p\n 

of / on the prime powers p k with k > 1 to obtain an Euler product expansion 

g M^Mg TT P 3 TT P 3 ~P TT P 3 ~ P ~ 1 

m,n m 2 n 2 (p(mn) p 3 — 1 p 3 — 1 -f-^ p 3 — 1 

p\n p\m,p\n p\mn 

S ran -r-r (1 — p)p 3 -pr p 3 — p 



tt rr 

m 3 n 3 <j)(mn) p 3 — p — 1 11 , p 3 — p — 1 

p|n ' p|m,pfn 



m 3 ?! 3 1 1 ^3 _ p _ l 11 p3_p_2 

p|n p\m,p\n 

4.3. Corollary. Define r(n) as in Theorem 3. T/ien we have S\ jn = r(n)S and 
_V^V^ Kj) _ f To r ( n )S if n is odd; 



S T-r — p 4 -r-r p 3 + p 2 



□ 



2|i 



Proof. The first equality is immediate by taking m = 1 in 4.2. 

If n is odd, we have S' 2 n = S2, n = If 4 divides n, the condition 2|i 

in the definition of S 2 n is superfluous as p(j) vanishes for A\j. It therefore equals 
S 2,n = S hn = r(n)S. □ 

5. Proof of Theorems 2 and 3. 

We now have everything at our disposal to prove Theorems 2 and 3. 

Proof of Theorem 2. We substitute the value of [F^j : Q] from Lemma 3.2 into 
the expression for 8(a, b) provided by (2.2) to obtain 



oo oo 



(5.i) «(M) = EE/y " {J) 



By 3.2, there are only finitely many values of fij that can occur, namely the divisors 
of At. By 3.3, the value of fij only depends on the greatest common divisors of i 
and ij with certain integers depending on a and b. It follows that the set of pairs 
for which fcj equals a given divisor of At can be characterized in terms of a 
finite number of divisibility criteria on i and ij. This enables us to write S(a,b) as 
an integral linear combination of our basic sums S mjn for suitable values of m and 
n. As each of these sums is a rational multiple of S by 4.2, we conclude that S(a, b) 
is itself a rational multiple of S. 

It is not at all clear from the preceding argument that the resulting value for 
5(a,b) will always be positive. From the expression 5(a,b) = ^^=i^( a 7^) as a 
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value of i for which Si(a,b) is non-zero. If a is not a square, we can take i = 1 as 
8i(a, 6) = 5(a) is then positive by Hooley's result. For arbitrary a, there can be 
many values of i with 8i(a,b) = 0. In fact, for many i one can construct a that 
satisfy [F* : (a)] ^ i for almost all p. A list of such values of i can be found in [5, 
(8.9)-(8.13)]. The smallest value that is not in the list is i = 24, and we will show 
that not only 624(a), but also S24 (a, b) is always positive. 

We are interested in the primes p for which [F p : (a)] equals 24 and [F p : (6)] 
is divisible by 24. Up to finitely many exceptions, these are the primes that split 
completely in the field E = Q(C24, v^> v^), but not in any of its extensions 
E n = Q(C24n, 24 \/a, 2 \/b) for n > 1. By the results of Lenstra [5, Theorem 4.1], the 
set of these primes has positive density (under GRH) unless there is an obstruction 
'at a finite level', i.e., an integer h such that every automorphism a of the extension 
E C Eh is trivial on E n ^ for some divisor n(a) > 1 of h. Thus, it suffices to 
show that for each squarefree integer h, there exists a G Gal(^/Q) satisfying the 
following 2 conditions: 

1. (T is the identity on E; 

2. if p is a prime dividing h and q is the largest power of p dividing 24h, then a is 
not the identity on Q(Cq)- 

In order to construct such an automorphism, we observe that the maximal subfield 
E ah C E that is abelian over Q has the property that Ga,\(E ah / Q) is an elementary 
abelian 2-group. Assume without loss of generality that 6 divides h, and let q be 
a prime power as in condition 2. Take a q to be any non-trivial automorphism of 
Q(Cq) that is a square in Gal(Q(C 9 )/Q). As q is not a divisor of 24, the group 
Gal(Q(Cq)/Q) = (Z/qZ)* is not of exponent 2, and such an element a q exists. 
Define <7q as the automorphism of Q(C24/i) with restrictions ct|q(£ 9 ) = cr q . Then 
o"o is a square in Gal(Q((24h)/Q), so it is the identity on E n Q(C24/i) C E ah . 
This implies that there is a unique extension of a® to E(C,2±h) that is the identity 
on E. Any extension a of this automorphism of E(C,24h) to Eh now meets our 
requirements. □ 

Proof of Theorem 3. In the case where a, b G Q* are multiplicatively independent 
and Q*/ (—1, a, b) is torsionfree, we can use the values of fij from Lemma 3.2 and 
rewrite (5.1) explicitly as an integral linear combination of the type encountered in 
the proof of Theorem 2. The sums of the type S' 2 n from 4.3, which single out the 
contribution to S\, n of the terms with even i, can be used to obtain the compact 
expression 

S(a, 6)61,1 + S 2i A(a) + ^2,A(6) + S 2,A(ab) + ( 5 'l,lcm(2, A(o)) ~ «S ! 2,lcm(2,A(o))) ' 

The sum of the three terms involving A(a) yields S'i,i cm (2,A(a))) since we have 
S'2 A(o) = ^2 icm(2 A(o))' ^ s i mme diate if A(a) is even, and if A(a) is odd the 
equality S' 2 A ^ = S' 2 2 A( a ) f onows directly from the definition of the sum S' 2jn in 4.3. 
The explicit value of c a ,b = 5(a,b)/S now follows easily from the two statements 
in 4.3. □ 

6. Reformulation in terms of recurrent sequences 
If we write the rational numbers a and b in Theorems 2 and 3 as a = a\/a2 and 
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sequence {6 2 a^— b\a^ }^L has positive density. This formulation in terms of integers 
is useful in proving the unconditional result in Theorem 1. The proof given below, 
which is entirely elementary, is an easy extension of the argument of Polya occurring 
in [11, Chapter 8, problem 107]. A generalization to integer sequences of the form 
{X)i=i c i a ?}n^=o f° r arbitrary k > 2 was given by Polya in [10]. 

Proof of Theorem 1. Write a = a±/a 2 and 6 = 61/62 as above, and take 
gcd(ai, 02) = gcd(6i, 62) = 1- We have a± 7^ ±02 by the hypothesis that a, b G Q* 
be multiplicativily independent, so \x n \ tends to infinity with n. We need to show 
that the set S of primes that divide x n = 62a" — 61 a% for some n > is infinite. 

Suppose that S is finite, and set I = <fi(\xo\ ■ Yl P esP)- Clearly, we have t > 0. 
We derive a contradiction by showing that the sequence {x£ n }^ =0 is bounded. As 
£ is finite, it suffices to show that ord p (xe n ) remains bounded as a function of n 
for each p G S. Suppose that p G S is a prime that does not divide a±a2- Then we 
have ord p (a;^ n ) = ord p (xo) for all n since 

x £n -x = b 2 (a{ n - 1) - b 1 (a e 2 n - 1) 

is by the definition of £ divisible by p ord p( x o)+i Suppose that p G S is a prime 
dividing a\a2- Then p divides exactly one of a± and a2, say a\, and we have 
ord p (x£ n ) = ord p (6i) for all sufficiently large n. □ 

Integer sequences of the form {6 2 a^ — ^la^j^Lo are nnea r recurrent sequences of 
order 2. They can be defined by the recursion Xk+2 = (&i + 0' 2 )xk+i — aia^fc for 
all k > and the initial values xq = 62 — b± and x± = 62^1 — b± <i2- 

Much effort has been spent on the determination of the set of primes dividing 
linear recurrent integer sequences, see [2] and the references given there. In the case 
of second order sequences, our Theorems 1-3 lead to the following result. 

Theorem 4. Let r, s G Q be rational numbers, and 1Z = {x n }^ =Q an integer 
sequence satisfying the second order recursion Xk+2 = rxk+i — sxk for all k > 0. 
Suppose that X 2 — rX + s splits in Q[X], and that TZ does not satisfy a first order 
recursion. Then the set of primes that divide some term of 1Z is infinite; if we 
assume the generalized Riemann hypothesis, it has positive density. 

Proof. Let 01,02 £ Q be the two roots of the polynomial X 2 — rX + s. 

Suppose first that we have ai 7^ 02- Then we have x n = 62a" — 6102 for certain 
61,62 G Q. As X does not satisfy a first order recurrence, we have ai<2 2 7^ and 
6i6 2 7^ 0. After replacing, if necessary, the sequence {x n }^ =0 by {X n x n }^ ) =0 for 
suitable A G Q*, we may assume that ai and 02 are coprime integers. This only 
changes the set of primes that divide some term of 1Z by finitely many primes. 
In a similar way, after replacing {x n }^ =0 by {Xx n }^ =0 for suitable A G Q*, we 
may assume that 61 and 62 are coprime integers. If a = a\/a 2 and 6 = 61/62 are 
multiplicativily independent in Q*, we are in the situation of Theorems 1 and 2, 
and we are done. If a and 6 are multiplicativily dependent, then the results on 
torsion sequences from [14] imply unconditionally that the set of primes that divide 
some term of 1Z has positive density. 

Suppose next that we are in the inseparable case a\ = a 2 . Then we have x n = 
(61 + 62^)0" for certain 61, 62 G Q, and 62 7^ by assumption. Now all primes that 
do not divide 62 divide some term of 1Z, so we obtain a set of prime divisors of 
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The hypothesis that 1Z does not satisfy a first order recursion in Theorem 4 is only 
there to exclude trivialities. In order to remove the assumption that X 2 — rX + s 
splits in Q[X], one needs to prove the analogues of our Theorems 1 and 2 for the set 
(1.2) in the case where a and b are elements of norm 1 in a quadratic number field 
K and F p is replaced by the ring of integers of K modulo the principal ideal (p). 
It turns out that the inert primes lead to various complications. The torsion case 
can be found in [14] and, for the special case where a comes from the fundamental 
unit of K, in [7]. For a treatment of the non-torsion example proposed by Lagarias 
[4, p. 451], we refer to [8]. 

7. Numerical data 

Just like Artin's constant A = prime (1 — l/p(p — 1)), the universal constant 
S = Yip primed - p/ (P 3 ~ 1)) i n Theorem 2 is defined by a slowly converging product. 
One can obtain good numerical approximations to S, such as the approximation 

S « 0.57595 99688 92945 43964 31633 75492 49669 25065 13967 17649 

up to 50 decimal digits, by expressing — log S as a rapidly converging series involving 
the values ((d) of the Riemann zeta-function at arguments d > 2. This is done for 
Artin's constant in [1], and for S the expression 



oo , 

log S = - log C(3) + E lo S C(d) ^ ( ' 

A — O U I A \ 



d=2 k\d 



is derived in [6]. Here is defined by its initial values a\ = 0, = 2, as = 3 and 
the recursion formula cik+3 = flfe+i + a k f° r k > 1. 

It is not computationally feasible to determine the rational numbers c a ^ in The- 
orem 2 from numerical data. In the torsionfree case occurring in Theorem 3, the 
value of c a ,b lies between the extremal values 

9343 , 28001 

c m in = c 2 ,5 = 77^77 « -981 and c max = c 5 , 3 = _ _ n w 1.023. 

For the 41535 primes contained in the interval [7,500 000], one finds that we have 
5 G (2 mod p) for 23498 primes and 3 G (5 mod p) for 24429 primes. When divided 
by S, these fractions are approximately equal to .9823 and 1.0212, respectively. This 
shows that the deviations from S are 'numerically visible' in a qualitative sense, 
but it also makes clear that one cannot determine the fraction c a ^ from such data. 
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